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FOREWORD 


This  is  the  final  report  on  work  done  by  United  Technologies  Research 
Center,  East  Hartford,  Connecticut,  for  the  Air  Force  Avionics  Laboratory 
(DHO),  Air  Force  Systems  Command  (AFSC),  under  Contract  F336l5-75-C-1047> 
Project  2001-01-35  during  the  period  January  1975  to  December  1975*  The 
work  was  monitored  for  the  Air  Force  by  Mr.  0.  P.  Breaux,  AFAL/DHO.  The 
report  describes  the  results  of  an  investigation  to  demonstrate  a  vibrational 
energy  transfer  laser  operating  in  the  3  to  5  micron  wavelength  band. 

The  author  wishes  to  acknowledge  the  assistance  of  Mr.  R.  J.  Hall  who 
performed  the  numerical  modeling;  Mr.  V.  Failla  who  assembled  and  operated 
the  experiments;  and  Mr.  W.  L.  Nighan  and  Dr.  C.  J.  Ultee  for  helpful 
discussions. 

The  authors  submitted  this  report  for  Air  Force  review  on  3  February 

1976. 
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SECTION  I 
SUMMARY 


dm* 


The  vibrational  energy  transfer  (VET)  lasers  have  been  known  for  over  ten 
years.  In  this  class  of  device,  pumping  is  based  upon  the  excitation  of  the 
vibrational  energy  levels  of  one  molecular  species,  typically  a  diatomic  gas, 
and  the  rapid  exchange  of  those  vibrational  quanta  to  a  second  selected  molec¬ 
ular  species,  typically  a  known  laser  emitter,  by  near  resonant  energy  trans¬ 
fer.  This  technique  has  only  recently  received  attention  for  applications  in 
the  mid- infrared,  which  would  fulfill  an  important  Air  Force/DoD  requirement. 

Computer  modeling  of  these  systems,  conducted  at  the  United  Technologies 
Research  Center,  using  reliable  physiochemical  data  and  an  accepted  analysis 
of  electron/molecular  and  molecule/molecule  interactions,  has  indicated  that  a 
cw  electric  discharge  excited  vibrational  energy  transfer  mixing  laser  can  be 
an  efficient  source  for  continuous  laser  radiation  in  the  desired  4.0to5.0|j,m 
emission  band. 

Experiments  have  been  conducted  on  a  laboratory  VET  laser  device  utilizing 
the  Do  */ HC1  system  which  has  been  identified  as  having  the  most  promising 
potential  in  the  mid- infrared.  In  these  tests,  spontaneous  emission  from 
transfer  excited  HCl  has  been  observed  and  is  reported  herein.  The  degree  of 
excitation  in  HCl  has  been  deduced  from  these  measurements.  In  addition  the 
degree  of  vibrational  excitation  within  hydrogen  issuing  from  a  dc  discharge 
has  been  directly  measured  by  Raman  scattering.  These  measurements  indicate 
that  first  three  vibration  levels  of  HCl  excited,  in  this  manner,  are  charac¬ 
terized  by  a  vibrational  temperature,  Ty  »  1400  -  196o°K.  The  post  discharge 
D2  vibrational  temperature  implied  by  these  determinations  indicated  T^  = 

1300  -  1700°K.  Hie  Raman  scattering  measurements  indicated  T^  =»  900  -  138o°K 
characterizing  the  state  of  hydrogen  downstream  of  a  similar  discharge.  Vibra¬ 
tional  relaxation  has  been  measured  in  these  experiments  and  is  much  higher 
than  anticipated. 

The  relatively  low  discharge  efficiencies  suggested  by  these  measurements 
have  been  confirmed  by  calorimetric  heat  balance  measurements  made  on  water 
cooled  discharges.  The  exact  explanation  of  the  low  discharge  efficiency  has 
yet  to  be  determined.  Titration  experiments  using  nitrosyl  chlorides  to  mea¬ 
sure  post  discharge  atom  concentrations  indicate  concentration  levels  below 
the  sensitivity  of  the  experiment,  as  predicted  by  theoretical  calculations. 
Therefore,  it  is  concluded  that  H-atom  deactivation  of  the  donor  gas  species 
is  not  of  sufficient  magnitude  to  account  solely  for  the  low  degree  of  excita¬ 
tion  observed. 
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SECTION  II 


INTRODUCTION 


Objective 

For  certain  important  Air  Force  and  DoD  missions,  an  efficient  low  to 
medium  power  (10  to  100  watt)  source  of  infrared  laser  radiation  is  required. 
Of  particular  interest  is  the  4.0  to  5.0  p,m  region  of  the  spectrum,  corre¬ 
sponding  to  an  atmospheric  transmission  window.  One  of  the  more  promising 
classes  of  devices  to  satisfy  this  requirement  is  the  cw  vibrational  energy 
transfer  (VET)  laser.  The  objective  of  the  research  reported  here  is  to 
develop  the  technology  base  sufficiently  to  demonstrate  the  technical  feasi¬ 
bility  of  the  VET  laser  and  determine  the  potential  performance  characteris¬ 
tics  of  such  devices  for  airborne  use. 


The  Vibrational  Energy  Transfer  Laser 

The  VET  laser  operates  essentially  in  two  sequential  steps:  (1)  a  doner 
gas,  preferably  a  diatomic  molecule,  is  vibrationally  excited  in  a  separate 
excitation  zone.  (2)  A  second  gas,  the  lasing  species,  is  intermixed  with  the 
donor  gas,  resulting  in  the  transfer  of  vibrational  quanta  from  the  excited 
donor  to  the  laser  species.  With  carefully  matched  donor/acceptor  molecule 
pairs  and  carefully  chosen  operating  conditions,  laser  emission  from  the 
acceptor  species  can  be  sustained  continuously  and  efficiently  via  this  pump¬ 
ing  route. 

There  are  several  distinct  advantages  of  the  VET  laser  concept.  The 
region  of  vibrational  excitation  can  be  separated  physically  from  the  region 
of  laser  action.  Thus  different  conditions  can  be  established  in  each  of 
these  different  regions,  affording  independent  optimization  of  the  excitation 
step  and  the  transfer  and  lasing  step.  A  laser  species,  which  either  could 
interfere  with  or  could  be  decomposed  in  the  excitation  process,  can  be  removed 
from  the  excitation  zone  of  the  VET  laser,  Therefore,  a  much  greater  range  of 
choice  of  laser  species,  and  hence,  a  much  greater  choice  of  laser  emission 
wavelengths ,  is  created.  With  a  single  donor  molecule,  the  laser  emission 
frequency  can  be  tuned  by  selection  of  different  acceptor  species,  any  of 
which  could  be  pumped  by  the  single  donor  molecules.  Certain  donor  molecules 
which  can  be  excited  efficiently,  but  do  not  emit  in  the  IR  (homonuclear 
diatomics)  can  be  utilized  very  effectively  in  the  VET  laser  device.  Vibra¬ 
tional  quanta  are  transferred  from  the  donor  species,  which  do  not  exhibit 
optical  gain,  to  an  acceptor  molecule  which  can  exhibit  optical  gain  (examples 


include  heteronuclear  diatomics,  triatomics ) .  If  the  vibrational  energy- 
levels  of  the  donor  and  laser  species  are  closely  matched,  the  transfer  of 
vibrational  quanta  between  these  two  species  can  be  very  fast  and  very  effi¬ 
cient.  Such  near-resonant  vibrational  energy  transfer  has  formed  the  basis 
for  the  pumping  of  the  well-known  Ng-CC^  laser. 

Based  on  previous  calculations  (Ref.  l)  a  dc  electric  discharge  appears 
to  be  the  most  promising  excitation  source.  Within  an  optimized  glow  dis¬ 
charge,  electrical  energy  can  be  coupled  very  efficiently  into  the  vibrational 
modes  of  diatomic  gases,  and  thus  allows  for  the  study  of  many  different  prom¬ 
ising  donor  gases.  This  rather  broad  range  of  study  is  reported  herein. 
Because  the  transfer  and  lasing  zones  are  physically  separate  from  the  elec¬ 
tric  excitation  in  the  proposed  VET  laser  device,  a  similarly  broad  range  of 
laser  species  can  be  explored. 

Calcinations  indicate  that  the  use  of  donor  (starred)/emitter  pairs  such 
as  D2*-HBr  and  D0*-HC1  in  the  cw  electric  discharge  VET  laser  can  yield  small- 
signal  gain  coefficients  >1.0  percent/cm  and  specific  laser  output  powers  & 

50  kj/lbm  on  a  series  of  vibrational-rotational  transitions  of  the  hydrogen 
halide  molecules  distributed  through  the  4.0 -5.0  (j,m  region.  These  same  com¬ 
puter  calculations  (Ref.  1)  have  identified  the  little  recognized  sensitivity 
of  the  VET  laser  medium  to  poisoning  by  small  concentrations  of  free  atoms. 
Concentrations  of  free  atoms  (e.g.,  H,  F,  Cl)  as  low  as  1.0  mol-percent  are 
expected  to  deplete  optical  gain  in  the  VET  laser  medium  completely  via  the 
extremely  rapid  relaxation  of  the  vibrationally  excited  laser  emitter,  envi¬ 
sioned  to  be  a  hydrogen  halide.  This  very  rapid  vibrational  energy  decay 
induced  by  atom  collisions  with  the  emitter  completely  dominates  over  vibra¬ 
tional  energy  transfer  pumping  of  the  emitter.  Fortunately,  continuous  self- 
sustained  dc  glow  discharges  can  be  maintained  in  selected  diatomic  gases 
(e.g.,  Hp,  D2)  with  negligible  free  atom  production. 

The  experimental  appartus  and  approach  chosen  for  the  current  investigation, 
is  designed  to  elucidate  the  vibrational  energy  transfer  mechanisms,  the  effi¬ 
ciency  of  donor  excitation  and  the  degree  of  excitation  in  the  laser  gas.  The 
experiment  thus  has  been  designed  to  yield  maximum  scientific  information. 

The  final  airborne  version  of  the  VET  laser,  however,  may  be  considerably  dif¬ 
ferent.  The  electric  discharge  VET  laser  may  be  adaptable  to  premixed,  closed- 
cycle  operation.  Efficient,  light-weight  electrical  power  supplies  matched  to 
the  power  requirements  envisioned  for  the  device  are  available  and  could  be 
easily  integrated  into  the  airborne  design.  A  fully  developed  electrically- 
pumped  VET  laser  could  be  envisioned  to  operate  in  a  closed-cycle,  thereby 
obviating  storage  and  supply  problems.  In  this  mode  of  operation,  an  initial 
charge  of  donor  and  acceptor  gases  would  be  pumped  continuously  through  a 
closed  loop  provided  with  the  necessary  discharge  and  cooling  section.  Such 
a  system  shows  great  promise  for  airborne  applications. 
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SECTION  III 

EXPERIMENTAL  INVESTIGATION 


The  United  'Technologies  Research  Center  for  the  past  year  has  investigated 
the  VET  laser  concept  as  applied  specifically  to  the  DgVHCl  system.  This 
investigation  has  had  as  an  objective  the  study  of  the  characteristics  of  an 
actual  device.  The  motivation  for  the  approach  is  that  by  characterizing  the 
conditions  pertaining  in  an  actual  device,  (l)  an  assessment  can  be  made  of 
the  appropriateness  of  theoretical  models  of  the  VET  laser  system  which  have 
been  developed,  (2)  the  laser  potential  can  be  predicted,  and  (3)  a  practical 
knowledge  can  be  obtained  about  some  of  the  major  technology  problem  areas  (at 
least  those  which  are  common  to  an  airborne  and  a  laboratory  environment) 
which  must  be  overcome  before  the  concept  can  be  reduced  to  practice.  The 
primary  diagnostic  tool  utilized  in  this  study  has  been  the  observation  of 
spontaneous  emission  from  vibrationally  excited  laser  species  formed  as  a 
result  of  mixing  with  the  discharge  excited  species.  In  addition,  direct 
measurements  have  been  made  of  the  vibrational  excitation  within  the  donor 
species  downstream  of  a  self  sustained  dc  discharge.  This  experiment  will  be 
described  in  this  section. 

Electric  Discharge  Tubes  for  Donor  Gas  Excitation 

The  self  sustained  electric  discharges  used  in  this  study  were  of  a  simple 
design.  The  discharge  tubes  are  15  cm  long  with  a  2.54  cm  ID.  The  central 
portion,  10  cm  in  extent,  is  cooled  by  means  of  a  coolant  circulation  through 
a  concentric  jacket.  The  coolant  is  either  water  or  liquid  nitrogen.  Elec¬ 
trode  geometry  used  in  the  majority  of  experiments  in  this  study  consisted  of 
an  upstream  electrode,  which  was  a  0.6  cm  tube  with  slightly  flared  end,  and  a 
downstream  electrode  (cooled)  which  was  a  disc  penetrated  by  a  2.0  cm  hole, 
with  rounded  edges.  The  donor  gas  entered  the  discharge  through  the  tube  of 
the  upstream  electrode,  the  outside  of  which  was  insulated,  and  excited  through 
the  hole  in  the  downstream  electrode.  Under  most  conditions  the  flow  was 
sonic  or  near-sonic  at  the  exit  of  this  tube.  The  transit  time  through  the 
discharge  was  approximately  3  milliseconds.  The  discharge  circuit  was  bal¬ 
lasted  (20  kQ)  to  provide  stability. 

The  discharges  were  operated  in  the  range  of  p  =  10  -  20  Torr,  with 
discharge  currents  from  30  to  150  mA  through  each  tube.  The  typical  Townsend 
parameter  electric  field  strength-to-gas  number  density,  was  E/n  a  2.0  x  10^  - 
2.5  x  lO1^  V-cm2.  At  low  pressure  the  presence  of  striations  in  these  blue 
discharges  was  pronounced.  Striations  were  not  apparent  to  the  eye  at  pressures 
typical  of  actual  operating  conditions.  No  signs  of  arcing  were  observed, 
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however,  the  liquid  nitrogen  cooled  discharge  tube  appeared  to  have  a  central 
core  of  more  intense  visible  emission. 


Measurement  of  Vibrational  Temperature  by  Raman  Scattering 

Direct  measurement  of  the  vibrational  excitation  of  the  most  interesting 
donor  species  by  infrared  emission  techniques  is  precluded  by  the  lack  of  a 
dipole  transition  moment  of  the  homonuclear  donor  species.  Of  course,  this 
same  fact  insures  that  vibrational  energy  once  established  decays  relatively 
slowly.  Vibrational  excitation  of  homonuclear  species  can  be  measured  by 
Raman  scattering  because  the  occurrence  of  the  Raman  spectrum  depends  on  the 
molecular  polarizability  and  not  on  the  presence  of  a  permanent  dipole  moment. 

In  the  Raman  effect  a  photon  incident  on  a  gas  molecule  is  inelastically 
scattered.  The  scattering  process  can  be  thought  of  as  an  absorption  and 
simultaneous  reemission  from  a  virtual  state.  Stokes  shifted  emission  is 
observed  when  the  incident  photon  excites  the  molecules  to  a  higher  vibrational 
state.  The  scattered  photon  energy  is  less  than  the  incident  photon  energy  by 
the  vibrational  quantum  energy  and  hence  the  outgoing  light  is  down-  or  red- 
shifted.  The  opposite  process  whereby  the  photon  gains  an  amount  of  energy 
corresponding  to  the  vibrational  quantum  energy  is  called  anti-Stokes  scatter¬ 
ing.  For  anti-Stokes  scattering  to  occur  the  molecule  must  be  in  a  vibration- 
ally  excited  state.  Under  usual  thermal  conditions  the  anti-Stokes  component 
is  not  observed  because  the  degree  of  vibrational  excitation  is  low  for  most 
gas  at  temperatures  less  than  1000°K.  However,  under  discharge  conditions 
where  the  electrons  are  coupled  more  effectively  to  the  vibrational  modes  than 
to  the  translational  modes,  the  vibrational  temperature  can  substantially 
exceed  the  translational  temperature  and  the  anti-Stokes  component  can  be  com¬ 
parable  or  larger  than  the  Stokes  component.  It  can  be  shown  that  the  ratio 
of  anti-Stokes  to  Stokes  scattered  intensities,  Ias  and  I  respectively  is 
given  by 


where  v0  is  the  incident  photon  frequency  (wavenumber)  and  vv  is  the  vibrational 
frequency  (wavenumber)  and  N(v=»l)  and  N(v  =  0)  are  the  populations  of  the  first 
vibrationally  excited  state  and  the  ground  vibrational  state  respectively. 

Thus,  the  degree  of  vibrational  excitation  induced  by  an  electric  discharge 
can  be  determined  by  measuring  the  ratio  of  anti-Stokes  and  Stokes  scattered 
light  intensities. 
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description  of  the  Raman  Apparatus 


In  principal  it  would  be  possible  to  determine  the  degree  of  vibrational 
excitation  of  the  donor  species  (and  the  acceptor  species  as  well)  directly  in 
the  transfer  apparatus.  However,  because  of  the  anticipated  difficulties  of 
doing  the  measurements  it  was  decided  to  conduct  the  Raman  scattering  tests  in 
a  separate  apparatus  where  collection  efficiency  and  stray  light  rejection 
could  be  optimized.  Figure  1  is  a  schematic  diagram  of  the  laser  Raman  scat¬ 
tering  apparatus .  The  beam  from  a  Spectra  Physics  (Model  166)  argon-ion  laser 
was  brought  to  a  focus  in  the  center  of  the  scattering  cell  by  lens  L^.  The 
divergent  beam  emerging  from  the  cell  was  absorbed  in  a  beam  dump  composed  of 
two  colored  glass  filters  set  at  Brewster*  s  angle  to  the  incident  beam  and 
inside  a  blackened  box.  Light  from  the  focal  volume  of  was  collected  at  90 
degrees  to  the  incident  laser  beam  by  an  achromatic  lens,  Lg  (f/l.7).  lens  Lg 
collimated  the  light  and  transmitted  it  to  lens  Lg,  being  turned  through  <~9Q 
degrees  at  a  plane  first  surface  mirror  which  was  used  to  facilitate  align¬ 
ment  .  Lens  Lo  refocussed  the  image  of  the  focal  volume  onto  the  entrance  slit 
of  a  Jarrell- Ash  (model  150)  1-m  monochromator  and  was.  chosen  to  match  the 
solid  angle  of  the  monochromator  entrance.  The  actual  orientation  of  the  laser 
beam  was  vertical  so  that  the  vertical  image  of  the  focal  volume  matched  the 
vertical  slit  of  the  monochromator.  A  spherical  front  surface  mirror  Mg  col¬ 
lected  back  scattered  light  and  focused  it  at  the  focal  volume  and  into  Lg. 

The  scattering  cell  consisted  of  a  rectangular  parallelepiped  of  black 
anodized  aluminum.  Three  mutually  perpendicular  hol^s  were  drilled  through 
the  center  of  the  cell.  The  laser  beam  entered  through  a  window  perpendicular 
to  the  cell.  The  bore  of  this  hole  was  threaded  to  reduce  scattered  light. 

The  windows  were  sealed  with  thick  0-riags  which  allowed  adjustment  of  the 
windows  to  achieve  parallelism  and  thereby  minimize  excess  laser  beam  scatter¬ 
ing  by  multiple  reflection  and  beam  walk-off.  The  second  hole  allowed  col¬ 
lection  of  scattered  light.  Discharged  gases  flowed  thorough  the  focal,  volume 
along  the  third  hole.  The  surfaces  of  these  holes  were  sand  bless ted  before 
the  cell  was  anodized  to  minimize  reflections. 

The  detector  system  consisted  of  a  EMI  6256s.  photomultiplier  tube  in  a 
Products  for  Research  (Model  IE- 114  B)  LNg  cooled  housing  at  the  exit  slit  of 
the  monochromator.  A  lens  imaged  the  exit  slit  onto  the  BIT  cathode.  A  SSR 
Instrument  photon  counter  (Jtodel  ll4o)  and  discriminator  were  employed  to 
process  the  signal. 

The  discharge  used  in  the  tests  was  identical  in  form  to  the  discharge 
used  in  the  transfer  apparatus.  The  only  difference  between  the  Raman  dis¬ 
charges  and  the  transfer  discharge  was  that  the  Raman,  discharge  tube  and 
cathode  were  water  cooled:  instead  of  LWg.  cooled.  A  Rayleigh  horn  was  employed 


in  a  bend  in  the  glassware  joining  the  discharge  tube  to  the  scattering  cell. 
The  distance  between  the  cathode  exit  and  the  laser  focal  volume  was  ~  19  cm. 

'The  overall  system  spectral  response  of  the  Raman  system  is  shewn  in 
Fig.  2.  Response  was  measured  with  a  tungsten-quartz-iodine  lamp  operated  at 
2400°K  as  determined  by  an  optical  pyrometer.  The  photon  counting  rate  as  a 
function  of  wavelength  was  normalized  by  the  blackbody  function  (Ref.  2)  at 
2400°K  corrected  for  dependence  of  emissivity  on  temperature  and  wavelength 
(Ref.  3).  The  shape  of  this  curve  is  determined  mainly  by  the  photomultiplier 
response.  The  Raman  scattered  light  is  polarized,  in  the  present  case  perpen¬ 
dicular  to  the  grating  grooves.  The  effect  of  a  polaroid  filter  on  the 
response  at  the  Raman  wavelengths  was  determined  and  is  also  shown  on  Fig.  2 
along  with  transmission  of  a  filter  used  for  some  measurements.  The  curve 
gives  the  response  to  natural  light. 

The  monochromator  counter  was  calibrated  with  a  Hg  lamp. 

Experimental  Results 

In  order  to  optimize  the  system  performance,  measurements  of  the  Raman 
Stokes  band  of  nitrogen  were  made  at  higher  than  normal  discharge  pressures. 
Usually  the  entrance  slit  and  exit  slit  width  were  set  equal.  The  entrance 
slit  width  was  determined  by  observing  the  effect  of  narrowing  the  slit  width 
on  the  Raman  signal.  The  slit  was  set  at  the  width  which  just  began  to 
decrease  the  transmitted  signal  when  the  slits  were  narrowed.  Normally  the 
slit  width  was  0.6  mm  corresponding  to  spectral  full  width  at  half  intensity 
of  0.5  nm.  This  width  was  somewhat  larger  than  the  image  of  the  focal  volume 
would  be  for  the  laser*  beamed  brought  to  a  focus  by  a  diffraction  limited  lens 
(collection  system  magnification  was  4.5). 

The  Raman  Stokes  shifted  signal  for  Ng  signal  is  shown  as  a  function  cell 
pressure  in  Fig.  3.  The  data  labeled  Stokes  are  measured  at  the  intensity 
peak  on  scanning  the  Raman  band.  The  background  signal  level  was  recorded  0.6 
nm  from  the  signal  peak,  and  the  dark  count  was  measured  with  the  monochromator 
entrance  slit  blocked.  Hie  dependence  of  dark  count  on  cell  pressure  is 
attributed  to  finite  time  required  by  the  EMT  to  recover  from  exposure  to 
lights.  Hie  background  signal  is  mainly  due  to  stray  Rayleigh  scattered  light 
with  a  small  component  of  cell  wall  scattered  light.  The  highest  signel 
recorded  at  760  Torr  was  16,000  counts  per  second.  It  was  found  that  when 
properly  adjusted  the  back  mirror  Mg  increased  the  signal  ~  60  percent. 

Figure  4  shews  the  hydrogen  sensitivity.  Note  that  the  laser  power  at 
488  nm  is  1.4  watts  in  this  case  compared  to  1.0  watts  for  the  nitrogen  data 
shown  in  Fig.  3.  The  lower  overall  sensitivity  to  hydrogen  is  mainly  attrib¬ 
utable  to  lower  photomultiplier  response  at  the  longer  Hg  wavelength  (the 
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Figure  4.  Hydrogen  Stokes  Sensitivity 
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hydrogen  Raman  cross  section  is  approximately  twice  the  nitrogen  cross  section, 
Ref.  4).  ^he  fact  that  the  background  signal  is  less  dependent  on  cell  pres¬ 
sure  probably  reflects  the  lower  Rayleigh  scattering  coefficient  in  hydrogen 
(Ref.  5). 

When  measurements  of  anti-Stokes  signals  were  attempted  under  discharge 
conditions  several  problems  arose.  The  first  problem  to  manifest  itself  in 
nitrogen  was  the  presence  of  a  Lyman  grating  ghost  at  a  wavelength  equal  to 
0.9  of  the  parent  488  nm  argon  ion  line.  This  places  the  ghost  ~  1  nm  from 
the  Ng  anti-Stokes  line.  This  line  actually  presents  no  problems  to  the  mea¬ 
surement  of  the  anti-Stokes  signal  (it  can  be  discriminated  against  by  a 
Schott  BG-1  filter)  except  on  scanning  it  can  be  mistaken  for  the  Raman  signal. 
The  second  problem  arose  because  of  the  large  background  signal  levels  present 
in  the  afterglow.  Background  and  anti-Stokes  signal  levels  were  such  that  the 
anti -Stokes  signal  was  not  apparent  on  scanning  the  monochromator  across  the 
spectral  region  of  the  band.  For  the  anti-Stokes  measurements  the  monochroma¬ 
tor  was  set  at  the  Raman  anti-Stokes  wavelength  and  a  signal  averaging  tech¬ 
nique  was  employed.  This  techniques  has  been  employed  previously  to  minimize 
background  effects  in  Np(v)  measurements  (Ref.  6).  The  Raman  signal  is  deter¬ 
mined  from  four  consecutive  counting  periods  during  which  the  counting  rate  is 
averaged  on  a  chart  recorder.  The  Raman  signal  is  obtained  from 


C  -  C-ul  +  CQ0  -  C1Q  -  CQ1, 


where  C-^  is  the  average  counting  rate  with  the  discharge  and  the  laser  on 
(and  hence  includes  the  Raman  signal  together  with  the  background  signal  from 
the  discharge  and  scattered  unshifted  light),  Cqq  is  the  counting  rate  with 
the  laser  and  the  discharge  off,  C^q  is  the  counting  rate  with  only  the  laser 
on  and  CQ1  is  the  counting  rate  with  only  the  discharge  on. 

It  should  be  noted  that  the  result  of  this  procedure  must  be  carefully 
interpreted.  For  the  anti-Stokes  signal  the  result  is  straightforward 
because  the  anti-Stokes  signal  is  positively  correlated  with  the  presence  of 
the  discharge.  In  the  case  of  the  Stokes  line,  however,  a  signal  decrement  is 
measured.  The  reason  for  this  is  that  Stokes  signal  is  present  with  the  dis 
charge  off.  Hence,  what  is  measured  by  the  above  averaging  technique  is  the 
effect  that  the  presence  of  the  discharge  has  on  the  Stokes  signal.  This 
effect  will  be  discussed  in  the  next  section. 
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The  Effect  of  Vibrational  and  Gas  Temperatures  on  Stokes  Scattered  Intensities 


Equation  (1)  gives  the  ratio  of  anti-Stokes  signal  to  the  Stokes  signal 
when  the  afterglow  region  is  characterized  by  a  vibrational  temperature  T  and 
a  gas  temperature  T.  This  equation  cannot  be  applied  directly  to  data  in 
which  the  Stokes  signal  has  been  measured  with  the  discharge  off.  The  data  in 
Figs.  3  and  4  must  be  corrected  for  the  effect  of  Ty  and  T  before  being  used. 

When  a  molecule  in  a  state  characterized  by  vibrational  and  rotational 
quantum  numbers  v  and  J  respectively  undergoes  a  transition  in  a  state  v* ,  J* 
the  intensity  of  Raman  scattered  light  can  be  shown  (Ref.  7)  to  be 


„  n  f.  E(v.O)  .  hqBnj( J+! ) 

v’,J'  Qv(Tv)Qr(Tr)  V  ^  L  kTy  kT 


(?) 


where  n  is  the  gas  number  density  and  gj  is  nuclear  spin  degeneracy,  and 
Or  are  the  vibrational  and  rotational  contributions  to  the  partition  function, 
b$,  is  the  Placzek- Teller  coefficient  (Ref.  7),  E(v,0)  is  the  energy  of  the 
v-th  vibrational  level  for  J  =  0,  B0  is  the  rotational  constant  of  the  ground 
vibrational  state  and  the  remaining  physical  constants  have  their  usual  mean¬ 
ings.  Those  parts  of  this  relationship  which  depend  on  the  gas  density  and 
molecular  parameters  can  be  factored  out  to  define  the  relative  line  intensity 


=  gl(2J+l)  hJ 


exp 


hcB0<J(j+l)/kT  J 


(3) 


Now  the  rotational  partition  function  can  be  approximated  in  the  temperature 
range  of  interest  by 

Qr  «  (2i+l)2kT 
2hcB0 

where  t  is  the  nuclear  spin.  So  that  equation  (3)  can  be  written 


J' 


gl(3J+1)9R 
(21+1 J' 


exp 


J(J+i)9r 

T 


(4) 
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where  0  R  =' hcBc/k  is  the  characteristic  rotational  temperature.  Restricting 
discussion  to  the  vibrational  Q  branch  (J1  =  J)  which  is  most  intense  the 
Flaczek-Teller  coefficient  can  be  written  (Ref.  7) 


J 

J' 


_ J(J+D  +  a 

(2J-l)(2J+3) 


(Q-branch) 


(5) 


where  aQ  =  (3-4pp)/4Pp  is  the  trace  scattering  coefficient  and  P^  is  the 
depolarization  ratio  for  polarized  light. 

The  effect  of  gas  temperature  on  the  relative  Q-branch  intensities  can 
be  interpreted  in  a  straightforward  manner  from  Equation  (4).  The  T~?  depen¬ 
dence  arises  from  the  combined  effect  of  gas  density  and  rotational  partition 
function.  The  exponential  gives  the  effect  of  temperature  on  the  individual 
rotational  lines.  The  Q-branch  transitions  in  nitrogen  are  more  closely 
spaced  in  wavelength  than  for  hydrogen  because  of  the  higher  moment  of  inertia. 
In  fact,  all  the  lines  of  Ng  Q-branch  fall  within  the  bandpass  of  the  mono¬ 
chromator.  Summing  over  all  lines,  the  rotational  partition  function  is 
recovered,  so  that  for  Ng  £  S°J  T”^.  For  hydrogen  the  various  Q  branch 
transitions  are  separated  Enough  that  some  lines  fall  outside  of  the  mono¬ 
chromator  bandpass.  Thus  increasing  temperature  has  the  tendency  of  shifting 
the  peak  of  the  line  intensities  out  of  the  monochromator  bandpass.  The  loca¬ 
tions  of  the  first  four  Hp  Q-branch  transitions  have  been  calculated  and  are 
shown  in  Table  1  for  a  488  ran  exciting  line. 


TABLE  1 

Hg  Q  BRANCH  WAVELENGTHS 


Transition 


Wavelength  (ran) 


v*=0. 


J  =  0 
J  «=1 
J  =  2 


J  =  3 


612.27 

612.05 

6ll.6o 

610.93 


Now  for  hydrogen  gj  =  3  for  the  odd  J  levels,  and  gj  **  1  for  the  even  levels 
and  the  nuclear  spin  is  I  ■  l/3,  see  Ref.  6,  The  depolarization  ratio  has 
been  measured  for  natural  light:  Pn  *  0.070-0.073  (Refs.  8  and  9).  The 
depolarization  ratio  for  polarized  light  is  related  to  the  natural  light 
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depolarization  by  0p  =  pn(2"pn^»  30  'tha't  ‘bhe  trace  scattering  coefficient  is 
a0  =  18.7-19*8.  Here  aQ  =  19  has  been  used.  The  relative  intensities  of 
several  Q  branch  transitions  have  been  calculated  as  a  function  of  gas  temper¬ 
ature  and  are  shown  in  Fig.  5.  Because  of  the  variation  of  nuclear  degeneracy 
and  the  exponential,  the  S^l  intensity  is  dominant.  A  simple  sum  of  all  five 
transitions  has  a  dependence. 

Now  the  Stokes  signal  under  discharge  conditions,  Cs(Td,  Tv)  can  be 
related  to  the  Stokes  signal  with  the  discharge  off,  Cg(T0): 


where 


and 


CS(D) 

f(T) 

f(T) 


Cs(0)f(T) 

QvC? 7T 


=  ^2  for  No 


and  use  has  been  made  of  the  fact  that  Qy(T  *=  TQ)  «  1.0.  TQ  is  the  initial 
gas  temperature  and  Tjj  is  the  post  discharge  temperature.  Equation  (1)  can 
now  be  written 


W  exp  ^ 


C„(To)f(T) 

CaB(TE)  Rs  \v0-vv  J 


(6) 


where  Ras  and  Rg  are  the  system  responsivity  at  the  anti-Stokes  and  Stokes 
wavelengths,  respectively  (see  Fig.  2). 

The  post  discharge  vibrational  temperature  can  be  calculated  from  Equation 
(6)  using  the  gas  temperature  factor  f(T),  or  directly  from  discharge  measure¬ 
ments  of  both  Stokes  and  anti-Stokes.  Both  of  these  procedures  have  been  used. 
It  was  found  that,  in  nitrogen  that  the  signal  averaging  procedure  resulted  in 
taking  the  differences  (small)  of  large  numbers,  due  to  the  presence  of  a 
large  background  signal  when  the  BG1  filter  wasn't  used.  This  led  to  a  large 
experimental  uncertainty  so  that  except  for  one  measurement  in  which  the  BG1 
filter  was  used,  the  nitrogen  data  has  been  reduced  using  Equation  (6)  along 
with  a  post  discharge  gas  temperature  calculated  from  known  heating  rates. 

The  background  signals  in  the  case  of  Hg  weren't  large  and  the  results  of 
using  Equation  (6)  with  extrapolated  discharge  off  data  and  using  discharge 
data  directly  were  in  good  agreement.  The  results  of  several  experiments  are 
shown  in  Table  2. 
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relative  intensity 


TABLE  2 


RAMAN  VIBRATIONAL  TEMPERATURE  MEASUREMENTS 


p 

id 

Vd 

Oas_-i 

cs  , 

Tv 

(Torr) 

(mA) 

(kV) 

(sec  ) 

(sec  ) 

0 

K 

10 

30 

2-50 

2b 

4o 

1710 

14 

30 

2.15 

42 

87 * 

1360 

14 

6o 

2.00 

51 

80* 

1558 

14 

75 

1.85 

70 

78* 

2290 

20 

30 

2.30 

4.5 

21.5 

926 

20 

6o 

2.15 

10.5 

20 

1080 

20 

TO 

2.00 

13-5 

19 

1142 

11.5 

30 

30 

19* 

1350 

20 

50 

3-10 

19 

36* 

1140 

21 

30 

3.90 

41 

29 

1310+ 

Mata  reduced  using  Equation  (6)  and  Figs.  3  and  4. 

+all  data  except  this  used  the  sonic  tube  anode  discharge  design.  This  exper 
iment  used  a  radial  injector  design  discussed  in  Section  IV. 


Vibrational  Energy  Transfer  Experiments 

The  vibrational  energy  transfer  apparatus  has  been  described  previously 
(Ref.  l)  so  that  here  only  the  salient  features  will  be  reiterated  and  modifi¬ 
cations  made  to  the  previous  apparatus  will  be  described.  In  the  transfer 
apparatus  (Fig.  6)  gases  excited  in  two  2.54  cm  diameter  discharge  tubes  flow 
through  constant  area  transition  sections  to  two  0.4  cm  high,  12.5  cm  wide 
channels.  A  secondary  gas  is  injected  into  the  main  flow  through  88-0.010 
inch  holes  drilled  in  a  single  row,  through  the  top  and  bottom,  channel  walls 
at  the  beginning  of  the  parallel  wall  central,  section,  see  Fig.  7.  The  gases 
exit  through  two  transition  sections  to  a  chemical  trap  charged  with  zeolite 
and  hence  to  a  large  (nominally  600  cfta)  vacuum  pump.  For  safety,  nitrogen  is 
added  to  the  flow  after  the  chemical  trap.  The  channels  are  joined  together 
at  the  center  and  the  sidewalls  are  closed  by  Brewster  angle  window  holders  in 
which  CaF2  windows  are  used.  The  windows  are  purged  with  an  inert  gas. 

The  modifications  which  have  been  made  to  the  apparatus  include  the 
addition  of  cooling  coils  to  the  outer  top  and  bottom  surface  of  the  nickel 
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channels.  Alcohol  is  circulated  through  these  coils  and  through  a  heat 
exchanger  coil  in  a  dry  ice/alcohol  mixture  by  two  centrifugal  pumps.  In 
another  modification  the  injector  manifold  has  been  constructed  of  copper  bar 
and  is  bonded  to  the  channels.  This  manifold  is  cooled  by  a  loop  in  the 
alcohol  coolant  circuit.  Indium  has  been  used  to  seal  the  discharge  tubes 
which  are  LNg  cooled.  With  these  seals  a  leak  rate  corresponding  to  0.01  per¬ 
cent  of  the  total  flow  has  been  maintained  over  a  period  of  many  thermal 
cycles  of  the  apparatus. 

HC1  Band  Emission  Experiments 

Radiation  from  the  HC1  vibrational  fundamental  band  has  been  examined  in 
the  transfer  apparatus  as  a  function  of  distance  downstream  frcm  the  sidewall 
injector.  A  LNg  cooled  lead  sulfide  detector  with  a  3  mm  x  10  mm  chip  was 
moved  parallel  to  the  flow  direction  on  a  stable  translation  stage.  A  narrow 
bandpass  filter  (3.6  +  0.6  u  )  in  front  of  the  detector  isolated  the  low-lying 
vibrational  transitions  of  HC1.  A  tuning  fork  chopper  modulated  the  sidelight 
emission  signal  prior  to  amplification  by  a  PARC  (Model  124)  tuned  amplifier. 

A  7  mm  diameter  aperture  4  cm  from  the  detector  chip  defined  a  narrow  (5°  half 
angle)  collection  angle. 

The  behavior  of  sidelight  emission  versus  position  was  examined  as  a 
function  of  several  parameters  including  deuterir'm  (and  hydrogen)  flow  rate, 
hydrogen  chloride  flow  rate,  and  discharge  power.  Typical  data  representing 
the  effect  of  HC1  flow  rate  are  shown  in  fig.  8.  In  these  data  the  deuterium 
flow  rate  was  held  constant.  By  plotting  the  logarithm  of  the  signal  inten¬ 
sity  against  distance,  an  exponential  decay  length  can  be  deduced  from  the 
falloff  in  signal  intensity  at  the  longer  flow  distances.  WLth  a  knowledge  of 
the  gas  flow  rates  temperature,  and  pressure  the  gas  density  and  velocity  has 
been  calculated  and  a  decay  time  has  been  deduced.  By  noting  the  positions  of 
the  signal  maxima  the  effective  transfer  time  can  be  calculated.  The  result 
of  this  data  reduction  is  shown  in  Fig.  9* 

These  data  can  be  interpreted  using  a  simple  model  of  laser  fluorescence 
involving  V-V  energy  transfer  and  VT  relaxation  developed  by  C.  B.  Moore 
(Ref.  10).  This  model  has  given  excellent  results  for  the  Ng-C0g  transfer 
system  (Ref.  11),  and  although  the  Dg-HCl  system  is  not  as  simple  (involving 
more  levels  and  larger  anharmonicities)  this  model  represents  logical  starting 
point  in  the  interpretation  of  the  present  data.  Consider  a  mixture  of  an 
initially  unexcited  gas  instantaneously  mixed  with  a  vibrationally  excited  gas 
which  is  allowed  to  relax.  If  the  coupling  is  strong  (fast  W  transfer  rate) 
the  vibrational  degrees  of  freedom  of  the  gases  will  rapidly  equilibrate  and 
the  two  species  will  return  together  to  equilibrium  with  the  surroundings. 
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The  period  of  mutual  equilibration  is  characterized  by  the  W  transfer  time 
and  the  period  of  relaxation  is  determined  by  the  overall  relaxation  rate.  In 
Moore's  model  analysis  the  individual  population  densities  n^  are  given  by 


.  Xit  „  XPt 

n  =  A  e  x  B.e  <=■ 

i  1  1 


(7) 


For  the  initially  unexcited  gas  Aj_  =  B^.  The  exponentials  and  Ag  are  both 
negative  and  are  given  exactly  in  Ref.  10.  When  the  W  transfer  time  is  much 
faster  than  the  relaxation  rate  then  exponential  factors  are  given  approxi¬ 
mately  by 

'’■I  ■  MV  +  V™1] 


and 
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vihere  the  rate  constants  are  associated  with  the  following  reactions 

Irp  * 

D2(l)  +  HCl ( 0 )  -  D  (0)  +  HCl(l)  +  AE 
ke 

HCl  ( 1 )  +  HCl  -»x  HCl  ( 0 )  +  HCl 

kip 

HC1(1)  +  E2  “>  HCl  (0)  +  HCl 

Jr 

D2(l)  +  HCl  D  (1)  +  HCl 
D2(i)  +  d2  -i2  d2(0)  *  d2 


The  time  to  the  peak  HCl  density  can  be  found  by  differentiating  Equation  (7) 
and  setting  the  result  equal  to  zero,  thus 


t 

m 


(9) 
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Equation  (8)  has  been  solved  for  rate  constant  values  taken  from  the  literature 
(Refs,  12-14) .  The  result  is  shown  as  the  solid  line  in  the  upper  portion  of 
Fig.  9«  Note  that  it  has  been  necessary  to  add  a  constant  =0.3  msec  to  the 
data  to  simulate  the  effects  of  mixing.  A  constant  is  reasonable  because  of 
the  nearly  constant  total  pressure.  An  estimate  can  be  made  of  the  laminar 
mixing  time  from 

t  — 

mix  tD 

where  1  is  the  mixing  scale  (set  equal  to  the  distance  between  injection  holes) 
and  D  is  the  binary  diffusion  coefficient  (15).  The  mixing  time  is  calculated 
to  be  0.6  msec  for  the  present  mixer.  Thus  the  mixing  time  required  to  pro¬ 
duce  agreement  with  the  experimental  data  is  not  unreasonable. 

The  rate  constant  ^  is  much  faster  than  so  that  considering  the 
behavior  of  HC1  population  the  exp  (*]t)  term  rapidly  vanishes  and  the  HCl 
emission  decays  as  exp  (Ag^)*  The  lower  portion  of  Fig.  9  can  be  interpreted 
as  a  representation  of  the  experimental  variation  of  -Ag  as  a  function  of  the 
hydrogen  chloride  partial  pressure.  Comparing  the  experimental  variation  of 
relaxation  rate  with  the  rate  calculated  from  Eq.  (8)  using  the  individual 
rate  constants  calculated  by  previous  investigators  (Refs.  12-13)  show  agree¬ 
ment  with  the  numerator  of  the  Ag  expression.  However,  the  denominator  is 
large  because  of  the  large  excess  of  deuterium  over  hydrogen  chloride.  Hence, 
the  overall  decay  rate  is  in  poor  agreement  with  the  decay  rate  of  the  simple 
model. 

Spectrally  Resolved  Sidelight  Emission 

Emission  from  the  vibration- rotation  transitions  of  HCl  was  observed  at  a 
position  approximately  2.0  cm  downstream  from  the  injector  corresponding  very 
nearly  to  the  location  of  peak  HCl  band  emission  (discussed  in  the  previous 
paragraphs).  The  UTRC  emission  spectroscopy  diagnostic  system  has  been 
described  in  detail  previously  (Ref.  l6).  The  system  is  built  around  a 
McPherson  0.3  meter  (Model  218)  scanning  monochromator.  Emission  is  imaged 
one-to-one  onto  the  entrance  slit  of  this  instrument  by  aim  focal  length 
toroidal  front  surface  mirror.  Radiation  exiting  the  monochromator  is  focused 
by  an  ellipsoidal  mirror  onto  a  InSb  (PV)  detector.  The  detector  signal  is 
preamplified  and  introduced  into  a  PARC  (Model  124)  lock-in  amplifier.  A  2000 
Hz  chopper  in  front  of  the  entrance  slit  modulates  the  radiation  and  provides 
a  reference  signal  for  the  phase  sensitive  amplifier.  Normally  the  entrance 
slits  were  250 n  wide  (and  6  mm  high)  corresponding  to  a  resolution  of  approx¬ 
imately  2  cm”1.  The  emission  apparatus  was  focused  through  the  CaFg  channel 
windows  to  a  point  in  the  center  of  the  channel  nearest  to  the  monochromator. 
The  top  and  bottom  channel  walls  limit  the  collection  aperture.  The  mono¬ 
chromator  system  was  calibrated  by  a  blackbody  source,  normally  operated  at 
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800°K,  set  near  the  apparatus.  When  a  calibration  was  required  an  aluminum 
front  surface  mirror  reflected  radiation  from  the  blackbody  into  the  emission 
apparatus.  The  response,  checked  on  each  run,  was  relatively  flat  over  the 
spectral  range  explored  here. 

The  emission  from  vibrationally  excited  hydrogen  chloride  has  been 
recorded  on  a  strip  chart  recorder  along  with  a  wavelength  reference.  A 
typical  spectrum  is  shown  in  Fig.  10.  The  vibration-rotation  transitions 
(labeled  Py(j)  =  v,  J-l  -  v-1,  J  and  Rv(j)  =  v,  J+l  -  v-1,  J)  of  the  v  =  1^  0, 
2  -*  1  and  3  2  bands  are  clearly  resolved.  The  R  branch  has  been  resolved 

but  is  not  shown.  The  intensities  of  the  v  =  3  “*■  2  transitions  are  near  the 
detection  limit  of  the  monochromator  system.  The  emission  signal  is  uncor¬ 
rected  for  the  system  spectral  response.  In  order  to  do  so  the  peak  signal 
voltage,  Vs,  is  normalized  by 


I 

q 


=  Ml 

VBB 


where  VBB  is  the  blackbody  voltage  response  of  the  system  at  wavelength  X  and 
NX  is  the  blackbody  function,  usually  expressed  in  w  cm”2  ^m”1  SR-1,  and  I  is 
the  reduced  signal  intensity. 

Analysis  of  Spectral  Data 

In  order  to  deduce  population  distributions  from  the  reduced  intensities 
it  is  necessary  to  know  the  radiative  transition  rates.  The  intensity  of  a 
transition  v* ,  J'  -*•  v,  J  is  given  by 


(10) 


where  AyJ  is  the  Einstein  transition  probability  (sec"-1-),  Nvtjt  is  the 
population  density  of  the  energy  level  characterized  by  quantum  numbers  v  and 
J  (em"3),  and  yJ  is  the  wavenumber  of  the  transition  (cm"1).  Now  if 

the  vibrational  mode  is  characterized  by  temperature  Tv  and  the  rotational 
modes  by  TR  the  population  of  level  v’J’  can  be  written 


N(2J'+1) 


exp 


Gn(v* )hc  By  *  he  J 1 ( J'+l ) 

kTy  ‘  kTR 


(11) 


where  and  0R  are  the  vibrational  and  rotational  partition  functions 
respectively,  N  is  the  total  gas  number  density,  GQ(v' )  is  the  energy  of  v' 
level  with  respect  to  the  vibrational  ground  state  (in  cm"1),  By  is  the 
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Figure  10.  Typical  D2  — HCI  Spontaneous  Emission  Spectra 


rotational  energy  spacing  of  the  v'  level,  and  (2J'+l)  gives  the  degeneracy  of 
the  «T'  level.  The  Einstein  transition  probability  can  be  written 


vM'  _  fiW'vtj'.vJ  |Rj'|ZF^jJ,s 
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vj 


J' J 


3h 


2J'h  1 


(12) 


where 


I  V 

|Rv 


is  the  pure  vibrational  dipole  transition  moment, 


v'  J' 

*Vj  is 


vibration-rotation  interaction  factor  and  Sjij  is  the  rotational  transition 
moment  (fjt  j  -■  J'  for  the  R -branch  and  S^.,j  =  J'-tl  for  the  P-b ranch ) . 


The  fundamental  band  vibrational  transition  dipole  moment  is  usually 
taken  as  proportional  to  v'  (Ref.  17),  i.e., 


v' 

'-1 


(13) 


The  vibration-rotation  interaction  factor  has  been  measured  from  absorption 
spectra  (Ref.  18) ,  the  result  is 

F(m)  =  1  +  .026  m 

where  in  =  J'  for  the  P-branch  and  m  =  -(J'+l)  for  the  P-branch. 

By  combining  these  expressions  and  rearranging  it  is  easy  to  show  that 
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J'(j'+l) 


(14) 


where  C  -  64m  c|Rq|  "N^I/q^,.  l/v  SFv'  is  the  normalized  intensity.  If  the 
analysis  is  restricted  to  the  P-branch:  J'  «  j-1  -*•  J,  Equation  (l4)  can  be 
written 
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In  C  = 
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J(J-1) 


Thus  the  familiar  result,  the  logarithm  of  the  normalized  intensity  is 

proportional  to  J(J-1)  is  obtained.  The  rotational  temperature  can  be  obtained 

by  plotting  the  normalized  signal  intensities  of  the  various  transitions  P  (j) 

♦  * 
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as  a  function  j(j-l).  The  rotational  temperature  can  be  obtained  from  the 
slope  of  the  line  fitted  to  the  experimental  data.  The  vibrational  population 
densities  are  proportional  to  the  intercept  of  the  fitted  lines  with  the  axis. 
The  vibrational  temperature  can  then  be  deduced  from 
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I  G0(v')hc 

4  - - 

v  jFv'  intercept  klIv 


A  typical  rotational  distribution  plot  is  illustrated  in  Fig.  11.  The 
rotational  temperature  from  these  data  is  TR  =  3?5°K.  Note  that  same  transi¬ 
tions  are  missing,  specifically  P^(9)  and  P_(6),  because  they  overlap  each 
other  and  are  totally  unresolved.  The  scatter  of  the  3  -*•  2  data  is  typical 
and  is  probably  due  to  the  uncertainties  associated  with  nhe  low  s/N  ratio  of 
the  data  for  this  transition. 


The  curvature  of  the  1  -*■  0  data  is  caused  by  self  absorption.  A 
computer  program  has  been  written  tc  correct  the  measured  intensities  for  self 
absorption.  Given  estimates  for  the  HC1  partial  pressure,  the  temperature  and 
the  optical  path  length  X,  the  computer  code  calculates  the  absorption  coeffi¬ 
cient  <*j  for  each  v  =  1  -*■  0  P-branch  transition.  The  measured  intensity  IE(J) 
is  then  corrected  for  optical  thickness  according  to  the  relationship  (Ref.  19) 
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The  computer  code  determines  the  rotational  temperature  from  these  corrected 
intensities  using  a  least  squares  fitting  technique.  The  rotational  temper¬ 
ature  is  used  to  recalculate  the  absorption  coefficients  and  hence  th8»- cor¬ 
rected  intensities.  The  procedure  is  repeated  until  the  corrected  rotational 
temperature  converges . 

A  typical  vibrational  distribution  for  the  same  data  as  Fig.  11  is  shown 
in  Fig.  12.  The  vibrational  temperature  deduced  from  the  low  lying  transi¬ 
tions  is  Tv  =  1700°K.  The  actual  vibrational  state  number  densities  can  be 
calculated  from  the  overall  system  response.  This  has  been  estimated  frcm  an 
average  of  all  the  spectral  data.  Figure  13  shows  the  vibrational  number  den¬ 
sities  obtained  in  this  manner.  The  Nc  level  is  calculated  from  the  HC1  par¬ 
tial  pressure. 

Several  experiments  have  been  conducted  under  various  operating  conditions, 
including  substitution  of  hydrogen  for  deuterium  and  the  use  of  He  diluent. 
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Figure  12.  Typical  HCI  Vibrational  Distribution 


The  emission  spectra  have  been  recorded  and  analyzed  according  to  the  above 
described  procedures.  The  results  are  given  in  Table  3. 

Measurements  to  Determine  H-Atcm  Concentration 


The  optical  specific  power  performance  of  the  D2-HCI  VET  laser  has 
been  predicted  (Ref.  1)  to  be  sensitive  to  D2  dissociation.  For  this  reason 
an  attempt  was  made  to  determine  H-atom  concentrations  exiting  from  the 
standard  design  electric  discharges  used  in  this  study.  The  principal 
effect  of  H  atoms  is  to  de-excite  HC1  vibrationally  excited  molecules 
(Ref.  20);  H  atoms  are  270  times  more  effective  collision  partners  than 
HC1  molecules  are  at  300°K.  However,  rapid  relaxation  of  ^(v)  has  also 
been  measured  (Ref.  21);  H  atoms  are  600  times  more  effective  than  H2. 

Several  techniques  have  been  used  to  quantitatively  determine  H-atom 
concentrations:  (1)  isothermal  calorimeter  (Ref.  22),  (2)  Wrede  gage 
(Ref.  23)  and  (3)  nitrosyl  chloride  (N0C1)  titration  (Refs.  24  and  25). 

The  last  mentioned  technique  was  chosen  for  these  experiments  because  it  is 
easier  to  interpret  the  results  and  it  is  capable  of  operation  in  the 
pressure  range  wherein  the  discharge  is  expected  to  be  maintained.  The 
reaction  mechanism  sequence  for  H  atoms  and  N0C1  is  as  follows  (Ref.  25): 


H  +  N0C1  —  NO  +  HC1  (1) 

H  +  NO  r-WM  -  HNO*  +  M  (2) 

HNO*  -  HNO  +  hv  (3) 

/ 

When  the  concentration  of  N0C1  is  much  less  than  the  H  atom  concentration, 

NO  molecules  generated  in  reaction  (l)  rapidly  react  with  excess  H  atoms  to 
give  excited  HNO  which  emits  radiation  in  several  bands  in  the  visible  red 
and  near  infrared  region  (Ref.  24).  When  the  N0C1  concentration  is  equal  to 
or  exceeds  the  H  atom  concentration  all  H  atoms  are  consumed  in  reaction  (l) 
and  none  are  available  to  react  with  NO  and  form  HNO*,  the  emission  stops, 
and  a  clear  titration  end  point  is  exhibited. 

* 

The  apparatus  used  in  these  experiments  is  shown  schematically  in 
Fig.  l4.  The  titration  experiments  weren't  carried  out  in  the  transfer 
apparatus  because  nitrosyl  chloride  is  reactive  and  corrosive  to  some  of 
the  materials  used  in  the  transfer  apparatus.  As  can  be  seen  from  the 
schematic  diagram,  hydrogen  gas  flows  through  a  standard  design  discharge, 
through  a  turn  in  a  light  trapping  elbow,  past  the  nitrosyl  chloride  mixer 
and  into  a  30  cm  long,  1.27  cm  i.d.  flow  tube.  Nitrosyl  chloride  and 


32 


.j 

.v  ? 


m 


HYDROGEN 

FROM 

FLOWMETER 


STANDARD 
WATER-COOLED 
DISCHARGE  TUBE 


TO  LN2  COOLED 
COLD  TRAP  AND 
VACUUM  PUMP 


LIGHT  TRAP 


CHOPPER 


NITROSYL  CHLORIDE 
FROM  FLOWMETER 


FILTER 


TO  SCOPE  AND 
TUNED  AMPLIPHIER 


Figure  14,  Hydrogen  Atom  Titration  Experiment 
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SUMMARY  OF  TRANSFER  EMISSION  EXPER 


condensible  reaction  products  are  trapped  in  a  liquid  nitrogen  cold  trap, 
and  after  an  experiment  are  neutralized  in  aqueous  NaOH.  Nitrosyl  chloride 
was  used  without  further  purification  from  cylinders  supplied  by  Matheson. 
The  flowrate  of  nitrosyl  chloride  is  measured  by  a  rotameter  type  flowmeter 
equipped  with  glass  and  tantalum  floats.  Nickel  plated  fittings  and  teflon 
tubing  have  been  used  throughout  the  N0C1  supply  system.  Emission  from 
HNO  was  detected  through  the  glass  flow  tube  wall  by  a  red  sensitive  photo¬ 
multiplier  (Hamamatsu  RU56 )  equipped  with  a  Wratten  88a  filter  to  block 
visible  emission.  The  emission  was  chopped  and  amplified  in  a  tuned,  phase 
sensitive  amplifier  (PAR  HR-8). 

The  results  have  been  completely  negative  with  no  detectable  HNO 
emission.  An  estimate  of  the  system  sensitivity  can  be  made  from  the 
intensity  of  band  emission  (Ref.  24),  the  sensitivity  of  the  PMT  (from 
manufacturer’s  specifications)  and  estimates  of  the  collection  solid  angle 
and  volume.  The  minimum  detectable  H-atom  concentration  can  be  shown  to 
be: 

4ttVs  1/2 

^min  “  2  ^hvI0N0VcncSRL^ 


where  Vs  is  the  minimum  detectable  signal  voltage,  hv  is  the  quantum  energy 
of  the  transition,  I0  is  band  intensity  (see  Ref.  24),  N0  is  Avogadro's  number, 
Vc  is  the  collection  volume,  Oc  is  the  collection  solid  angle,  S  is  the  PMT 
anode  sensitivity  in  A/w  and  RL  is  the  PMT  load  resistor.  With  Vs  =  mV, 

Vc  =  1  cm^  nc/4TT  =  .0025, Rl  =  10^  ohms, and  S  =  4  x  10^  A/w,  the  minimum 
detectable  hydrogen  atom  concentration  is  ~  8  x  10"11  mole/cm^.  This 
corresponds  to  .03  percent  of  the  typical  Hg  concentration. 

Several  additional  experiments  were  conducted  in  an  attempt  to  observe 
emission  from  HNO*.  The  discharge  tube  was  operated  deliverately  far  off 
the  experimental  design  point  in  order  to  attempt  to  maximize  H-atom  produc¬ 
tion.  Typical  pressure  in  the  discharge  tube,  under  these  conditions,  was 
3  Torr  and  the  hydrogen  flow  was  reduced  to  1,6  mi Hi moles /sec  corresponding 
to  a  discharge  transit  time  of  approximately  6  msec.  The  N0C1  flowrate  was 
varied  over  a  range  of  6  to  210  micromoles/sec  corresponding  to  from  0.3  to 
11.5  percent  of  the  hydrogen  flowrate.  Tests  were  conducted  in  which  the 
entering  hydrogen  stream  was  allowed  to  bubble  through  water  in  order  to 
saturate  the  hydrogen  stream.  The  walls  of  the  discharge  tube  were  carefully 
cleaned  and  coated  with  a  sulphuric  acid.  The  dc  discharge  was  replaced  by 
a  microwave  cavity  (Ref.  26)  enclosing  a  l/2  inch  o.d.  quartz  tube  which 
extended  down  to  the  mixer  when  the  light  trapping  elbow  was  removed.  These 
tests  and  techniques  were  tried  in  succession  and  in  combination  with  no  HNO 
emission  observed. 
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Several  possible  explanations  of  the  negative  result  have  been  considered. 
One  possible  explanation,  of  course,  is  that  the  H  atom  concentration  emanating 
from  the  dc  discharges  is  less  than  the  detection  limits  of  the  apparatus. 

N0C1  decomposes  thermally  (Ref.  27)  and  photochemically  (Ref.  28);  however,  NO 
is  a  decomposition  product  and  it  should  give  HNO*  emission.  The  nitrosyl 
chloride  might  be  dissociated  by  back  diffusion  into  the  discharge,  but  again 
HNO*  should  still  be  seen.  A  reaction  may  take  place  between  N0C1  and  foreign 
substances  in  the  tube.  Finally,  quenching  of  the  HNO  may  occur. 

If  the  dissociation  fraction  were  less  than  0.03  percent  it  should  have 
a  small  effect  on  laser  performance. 


Discharge  Heat  Balance  Measurements 


The  low  vibrational  temperatures  within  the  donor  gas  as  determined  by 
both  the  Raman  scattering  experiment  and  the  transfer /mixing  experiments 
implies  that  the  electric  discharges  are  much  less  efficient  than  anticipated 
(Ref.  1).  For  this  reason,  heat  balance  measurements  of  the  discharges  were 
undertaken.  It  was  decided  to  conduct  these  experiments  in  an  apparatus 
separate  from  the  Raman  and  transfer  experiments  in  order  to  facilitate 
measuring  all  heat  losses.  Furthermore,  in  order  to  simplify  the  measurements 
it  was  decided  to  water  cool  the  discharge  tube.  The  discharge  tube  is 
identical  to  standard  discharge  previously  used.  The  flow  tube  (length  =  38  cm) 
immediately  downstream  of  the  discharge  tube  expands  to  a  larger  diameter 
D  <=»  5  cm.  Heat  losses  to  the  discharge  wall,  cathode,  and  this  flow  tube 
are  individually  determined  using  water  calorimeters.  The  temperature  rise 
across  each  component  is  sensed  by  thermocouples  inserted  in  the  cooling  water 
circulating  to  each  component .  The  water  flow  rate  is  measured  by  a 
calibrated  flowmeter.  An  injector  immediately  downstream  of  the  flow  tube 
section  can  be  used  to  inject  a  secondary  gas  into  the  primary  stream.  A 
thermocouple  junction  is  inserted  in  the  flow  11  cm  downstream  of  the  injector. 
Several  styles  of  thermocouples  have  been  used,  including  a  bare  thermocouple, 
a  teflon  coated  thermocouple,  and  a  flow  stagnating  thermocouple.  No 
meaningful  difference  in  sensible  temperature  have  been  noted  with  the 
different  thermocouple  styles.  Estimates  (Ref.  29)  of  the  average  number  of 
collisions  that  a  single  molecule  being  swept  past  the  probe  will  make  with 
the  probe  are  much  less  than  usual  number  of  collisions  required  for 
accomodation  of  vibrational  energy  at  metallic  or  teflon  surfaces  (Ref.  30). 
Therefore, the  temperature  measured  by  this  thermocouple  probe  has  been 
interpreted  as  representing  the  gas  translational  temperature. 

Heat  balance  measurements  have  been  made  in  nitrogen  and  hydrogen. 

Using  nitrogen  in  the  discharge  it  is  possible  to  add  CO2  through  the  injector, 
and  measure  the  total  power  added  to  the  gas.  Self  relaxation  in  CO2  is 
fast  enough  to  thermalize  the  mixture.  Thus  the  total  power  convected 
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downstream  of  the  flow  tube  can  be  determined  from  the  gas  translational 
temperature  measured  by  the  thermocouple  probe.  The  results  are  interesting. 
In  the  case  of  nitrogen  50-60  percent  of  the  discharge  power  is  dissipated 
through  the  walls  (including  the  cathode  which  accounts  for  >0  percent  of  the 
wall  losses),  10  percent  goes  into  translation  and  30-U0  percent  is  in 
higher  energy  states  (predominantly  vibration).  The  overall  power  measured 
totals  approximately  95  percent  of  the  input  discharge  power.  The  30-40 
percent  of  the  discharge  power  attributed  to  higher  energy  states  is  the 
power  which  is  measured  upon  injection  of  the  CO2.  At  the  gas  flowrates 
used,  this  "vibrational  excitation  efficiency"  corresponds  to  a  power 
loading  of  4.7  to  6.3  kj/mole  which  in  a  Boltzmann  vibrational  energy 
distribution  corresponds  to  Tv  =  l600-l820°K.  These  results  agree  rea¬ 
sonably  well  with  the  Raman  experiments. 

The  results  of  the  H2  discharge  heat  balance  tests  are  less  conclusive. 
Contrary  to  the  nitrogen  discharge  results,  nearly  all  of  the  input  discharge 
power  is  measured  as  sensible  heat  by  the  thermocouples.  Typical  values  of 
the  fractional  discharge  power  dissipation  in  various  regions  of  the  flow 
are  as  follows:  discharge  tube  walls,  30  percent;  cathode,  20  percent; 
flow  tube,  20  percent;  and  30  percent  is  sensed  as  power  converted  down¬ 
stream  of  the  flow  tube.  It  appears  then  that  70  percent  of  the  discharge 
power  is  dissipated  through  the  cooled  walls,  and  most  of  the  remaining 
power  is  sensed,  presumably  in  translation,  downstream  of  the  discharge. 
Deactivation  of  viorational  energy  at  the  thermocouple  probe  appears  to  be 
negligible  unless  the  accomodation  coefficient  for  H2  is  markedly  higher  than 
that  for  N2.  The  fraction  of  the  discharge  power  in  vibrational  energy 
modes  downstream  of  the  discharge  implied  by  these  measurements  is  small;  in 
fact  less  than  the  uncertainties  (5  percent)  in  the  measurements.  The 
electric  discharge  coupling  efficiency  is  definitely  lower  in  hydrogen  than 
it  is  in  nitrogen  for  the  conditions  explored. 

Laser  Tests 


Attempts  have  been  made  to  observe  stimulated  eramission  from  D^-HCl  in 
the  transfer  apparatus.  The  tests  were  not  successful.  For  these  tests  the 
resonator  was  comprised  of  external  mirrors  with  Brewster  angle  CaF2  windows. 
The  mirrors  used  were  a  silver  enhanced  silicon  total  reflecting  mirror  and 
a  ~  98  percent  reflecting  partial  transmitter.  A  liquid  nitrogen  cooled  lead 
sulfide  cell  equipped  with  bandpass  filters  was  used  to  detect  the  presence 
of  stimulated  emission.  Spontaneous  emission  w  observed  as  confirmed  by 
blocking  the  mirror  farthest  from  the  detector,  vtfide  ranges  of  flow  and 
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discharge  conditions  were  explored  (not  all  independently)  as  follows t 


Dg  flow 
HC1  flow 
static  pressure 
cavity  velocity 
discharge  current 
discharge  voltage 


.013  *  *038  moles/sec 
0  -  *001  moles/sec 

3.5-5  and  6-10  Torr 
lCr  cm/sec  and  5  x  10^  cm/sec 
1500-400  Ma  (both  tube) 

2. 3-4. 8  kV 


Tests  were  performed  with  and  without  zeolite  material  in  the  absorbent  trap 
as  is  indicated  by  the  two  pressure  ranges  and  two  cavity  velocities.  Removal 
of  zeolite  afforded  higher  pumping  capacities. 


Spontaneous  emission  was  not  spectrally  resolved  during  these  tests  so 
exact  information  is  not  available  as  to  vibrational  state  populations; 
however,  it  is  suspected  that  negative  results  are  caused  by  low  inversion 
densities. 
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SECTION  IV 


CONCEPTUAL  DESIGN 


A  preliminary  screening  has  been  conducted  of  vibrational  energy 
transfer  (VET)  candidate  molecular  pairs  in  order  to  determine  their 
potential  for  airborne  infrared  laser  system  applications.  The  D2-HCI 
system  has  been  selected  as  most  promising  based  on  atmospheric  transmission 
characteristics,  cycle  configuration,  system  potential  size  and  weight, 
and  aircraft  compatibility.  The  electric  discharge  VET  approach  appears 
to  offer  higher  potential  for  airborne  applications  than  thermal  or 
chemical  VET  excitation.  Preliminary  VET  laser  performance  calculations 
indicate  that  maximum  performance  may  require  a  low  gas  temperature 
within  the  active  medium.  Previous  investigation  (Ref.  31)  has  shewn  that 
the  most  effective  means  of  cooling  subsonic  flow  laser  mixtures  is  by 
a  liquid  N2  heat  exchanger.  This  technique  is  being  explored  in  the 
present  investigation.  Reconfirmation  of  these  estimates  is  required  when 
additional  analytical  and/or  experimental  data  become  available. 
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SECTION  V 


DISCUSSION  OF  RESULTS  AND  CONCLUSIONS 


Measurements  have  been  made  of  the  vibrational  excitation  of  hydrogen 
downstream  of  an  electric  discharge.  Measurements  have  also  been  made  of 
the  vibrational  excitation  of  hydrogen  chloride  as  a  result  of  mixing  with 
discharge  excited  deuterium  (and  hydrogen).  A  relationship  between  the 
measured  HC1  vibrational  temperature  and  the  corresponding  D^  vibrational 
temperature  can  be  derived  from  a  consideration  of  the  W  transfer  reaction 

D2(v=1)  +  HCl(v--O)  ^  D2(v=0)  +  HCl(v=l )  +  AE. 


The  rate  constants  are  the  same  as  defined  previously  and  AE  is  the  energy 
defect.  As  discussed  above  the  mixing  D2  and  HC1  streams,  in  the  transfer 
apparatus,  rapidly  reach  a  mutual  vibrational  equilibrium  in  a  time  compa¬ 
rable  to  the  W  transfer  time.  At  equilibrium,  the  vibrational  temperature 
of  deuterium,  TV2,  is  related  to  the  HC1  vibrational  temperature,  Ty^  by 
detailed  balance.  Using  Boltzmann  factors  to  describe  the  vibrational 
excitation  it  can  be  shown  that 


6  AT, 

2  AE 

T  '  kT  ’ 


where  0j  and  02  are  the  characteristic  vibrational  temperatures  (0  =  G(l) 
hc/k)  of  HC1  and  D2,  respectively.  Since  AE/k  =  02-0]_,  this  expression  can 
be  written 
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This  equation  shows  that  at  equilibrium  the  vibrational  temperature  of  D2 
is  less  than  HC1  temperature. 

An  estimate  of  the  initial  deuterium  vibrational  temperature,  i.e., 
prior  to  mixing,  can  be  made  from  the  equilibrium  vibrational  temperature 
by  assuming  that  ther  is  no  loss  of  vibrational  quanta  during  W  equilib¬ 
ration,,  This  Is  a  valid  assumption  if  VY  energy  transfer  is  much  faster 
than  VT  relaxation,  as  is  usually  the  case.  The  conservation  of  vibra¬ 
tional  quanta  can  be  written 
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where  G2(v)  and  Gp(v)  are  vibrational  energies  of  level  v  referred  to  the 
ground  vibrational  state  of  I>2  and  HC1  respectively,  and  N2  and  N-j_  are  the 
number  densities  of  and  HC1  respectively.  The  initial  deuterium  vibra¬ 
tional  temperatures,  TVo,  deduced  in  this  manner  are  shown  in  Table  3* 


The  power  convected  out  of  the  discharge  in  the  vibrationally  hot 
deuterium  flow  was  calculated  from  TVq  and  the  known  flowrates.  The  dis¬ 
charge  power  input  was  calculated  from  the  discharge  voltage  (across  the 
electrodes)  multiplied  by  the  discharge  current.  Note  that  in  run  number 
5  only  one  discharge  tube,  the  one  closest  to  the  spectrometer,  was 
activated.  The  numbers  shown  are  applicable  to  the  two  discharge  tubes 
activated  and,  therefore,  are  comparable  to  the  other  data  in  this  table. 


The  data  of  Table  3  tend  to  indicate  that  it  may  be  that  the  discharge 
is  inefficiently  coupled  to  the  gas  vibrational  modes.  If  a  correction  is 
made  for  the  estimated  power  going  into  the  cathode  fall,  the  estimated 
fractional  power  transfer  to  vibration  is  20-25  percent.  These  results 
are  somewhat  questionable  due  to  the  indication  of  rapid  decay  of  HC1  vibra¬ 
tional  excitation.  However,  it  should  be  noted  that  the  Raman  experiments 
support  the  transfer  experimental  results,  in  that  a  low  degree  of  vibra¬ 
tional  excitation  is  indicated.  The  transfer  experiments  indicate  TVq  = 
1300-1700°K  while  the  Raman  scattering  tests  indicate  TVq  =  1000-1350  °K. 

The  low  measured  donor  gas  excitation  can't  be  soley  attributed  to 
inefficient  discharge  coupling  at  this  time.  Another  possible  explanation 
is  that  the  vibrational  excitation  is  adequate,  however,  vibration  deactiva¬ 
tion  rapidly  thermalizes  the  excitation  in  the  time  that  it  takes  for  the 
discharge  donor  gas  to  flow  from  the  discharge  to  the  region  of  excitation 
measurement.  If  this  latter  explanation  is  the  dominant  effect  then  since 
the  Raman  experiments  show  low  vibrational  temperatures  rapid  vibrational 
excitation  of  the  deuteri  urn  itself  must  be  considered. 


The  results  of  several  Raman  scattering  experiments  are  shown  in  Table 
4.  The  efficiency  of  discharge  in  converting  electrical  energy  to  vibrational 
energy  in  the  gas  his  been  calculated  from  the  Raman  determinations  and  the 
electrical  power  input.  The  vibrational  energy  per  mole  is  calculated 
assuming  that  the  measured  vibrational  temperature  is  valid  for  all  levels 
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TABLE  4 

SUMMARY  OF  FAMAN  EXPERIMENTS 
WITH  ELECTRIC  DISCHARGE  EXCITED  GASES 


Nitrogen 

l4  Torr 

Sonic 

Tube 

n  =  .013  mole/sec 

/cl 

(mA) 

Vd 

(v) 

PDIS 

(v) 

Tv 

(°K) 

Ev(Tv) 

(kj/mole) 

n 

(*) 

30 

2150 

64.5 

1360 

2.8 

56 

Go 

2000 

120 

1558 

4.2 

44 

75 

1850 

139 

2000 

8.0 

74 

Hydrogen 

20  Torr 

Sonic 

Tube 

n  =  .030  raole/sec 

(n$) 

vd 

(V) 

(SF 

T 

V 

(°K) 

Ev(Tv) 

(kj/mole) 

n 

W 

3° 

2300 

69 

926 

0.08 

3.5 

60 

2150 

129 

1080 

0.20 

4.7 

90 

2000 

180 

1142 

0.27 

4.5 

Hydrogen 

21  Torr 

Radial  Injector 

n  =  .030  mole/sec 

i 

V, 

Ptitg 

T 

Ev(Ty) 

R 

(m£) 

(v) 

(SF 

(°K) 

(kj/mole) 

w 

30 

3900 

117 

1200 

0.35 

8.9 

30 

3900 

117 

1310 

0.54 

13.8 

up  the  vibrational  ladder.  It  is  noted  that  the  nitrogen  measurements, 
at  least  at  the  higher  discharge  current,  are  quite  reasonable  in  com¬ 
parison  to  theoretical  predictions.  The  efficiency  in  hydrogen  is  quite 
low. 


These  results  can  be  compared  with  other  previous  investigations. 

There  are  only  a  limited  number  of  direct  measurements  of  vibrational  exci¬ 
tation  from  electric  discharges.  Nelson  et  al  (Ref.  32)  used  Raman  scattering 
to  determine  vibrational  excitation  in  N2  1 5  cm  downs  tree"  of  a  transformer 
excited  discharge.  They  found  a  vibrational  temperature  Tv  =  1950°K.  In 
later  measurements  on  a  similar  N2  discharge  in  the  pressure  range  p  =  100- 
150  Barrett  and  Harvey,  using  a  Fabry-Perot  interferometric  Raman  technique, 
were  able  to  determine  Tv  =  1930-2356°K  (Ref.  33).  The  only  published 
previous  measurement  of  vibrational  excitation  in  a  hydrogen  discharge 
utilized  vacuum  UV  absorption  of  vibrational  states  (Ref.  3*0  at  3  Torr 
pressure  excited  in  a  microwave  discharge.  The  workers  in  this  study  esti¬ 
mated  that  1-4  percent  of  the  hydrogen  was  vibrationally  excited.  This  cor¬ 
responds  to  Tv  =  1300-1900°K. 

The  Raman  scattering  measurements  of  vibrational  excitation  reported 
herein  are  in  agreement  with  the  broad  range  of  previous  determinations.  In 
both  N2  and  Hg  they  are  on  the  low  side  of  these  previous  measurements.  This 
may  be  indicative  of  some  particular  inefficiency  associated  with  the  present 
discharges . 

The  heat  balance  measurements  discussed  in  the  previous  section  bear 
out  the  conclusion  that  the  present  discharges  are  characterized  by  low 
vibrational  excitation  efficiencies.  The  discharge  efficiencies  estimated 
from  the  Raman  experiments  are  In  good  agreement  with  the  heat  balance 
measurements.  The  very  low  discharge  efficiency  indicated  by  the  Raman 
measurements  in  hydrogen  is  within  the  experimental  errors  that  are 
expected  irom  the  heat  balance  measurements.  It  is  therefore  not  surprising 
that  negligible  efficiency  is  implied  from  those  experiments. 

No  definite  conclusions  can  be  drawn  as  to  the  reasons  for  the  low 
discharge  efficiencies.  It  appears  at  the  present  time  that  deactivation 
by  atomic  hydrogen  doesn't  account  for  the  low  efficiencies.  The  modeling 
calculations  indicate  that  although  some  deactivation  occurs,  even  at  the 
low  atom  concentrations  predicted  (.0016  percent)  the  decrease  in  vibra¬ 
tional  temperature  (less  than  10  percent)  is  not  nearly  enough  to  give  the 
low  efficiencies  observed,  unless  the  efficiency  was  low  to  begin  with  at 
the  exit  of  the  discharge.  The  nitrosyl  chloride  titration  experiments 
indicate  that  the  atomic  concentrations  issuing  from  the  standard  discharges 
is  less  than  the  sensitivity  of  the  experiment  (  —  .03  percent).  Although 
the  tiration  experiment  is  not  sensitive  enough  to  confirm  the  predicted 
concentrations,  it  should  be  sensitive  enough  to  measure  atomic  concentration 
levels  which  would  be  deleterious  to  laser  power  extraction. 
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The  higher  heat  losses  to  the  wall  measured  with  hydrogen  compared  to 
nitrogen  may  reflect  the  higher  thermal  diffusive ty.  Some  inefficiencies 
are  naturally  expected  because  of  the  strong  cooling.  Additional  ineffici¬ 
ency  is  expected  because  of  the  short  discharge  tubes  (to  minimize  H-atom 
production)  and  the  relatively  higher  ratio  of  cathode  potential  fall  com¬ 
pared  to  the  total  voltage. 

The  present  study  has  answered  several  important  questions  concerning 
the  generation  of  stimulated  emission  by  vibrational  energy  transfer  in  the 
D2-HCI  system.  Key  questions  remain  concerning  the  exact  nature  and  cause 
of  the  low  discharge  efficiency  in  producing  vibrational  excited  molecules. 
Whether  the  inefficiency  can  be  overcome  by  brute  force,  i.e.,  increased  power 
loading,  or  more  sophisticated  techniques,  has  yet  to  be  determined. 

Increased  testing  and  combination  of  tests  such  as  the  Raman  scattering 
experiment  and  heat  balance  measurements  will  answer  these  questions. 

Numerical  calculations  have  been  made  to  simulate  the  Raman  scattering 
measurements  of  post  discharge  hydrogen  vibrational  excitation.  Particular 
attention  was  paid  to  the  evolution  of  the  hydrogen  vibrational  levels 
through  the  12  cm  discharge  region  and  the  19  cm  downstream  of  the  discharge 
exit.  This  corresponds  to  the  position  where  the  Argon  ion  laser  beam  is 
focused  into  the  hydrogen  afterglow. 

Calculations  made  previously  of  the  electron  energy  distribution  in  a 
uniform  discharge  (Ref.  1),  give  the  electron  excitation  and  dissociation 
rates,  average  electron  energy  and  drift  velocity  as  functions  of  the 
Townsend  parameter.  These  rates  are  used  along  with  models  of  vibrational- 
vibration  energy  transfer  and  vibrational- translation  relaxation  by  molecular 
and  atomic  hydrogen.  The  buildup  of  Hp(v)  is  calculated  in  this  manner 
through  the  discharge.  At  the  discharge  exit  the  electron  collision  rates 
are  set  equal  to  zero  and  calculation  is  carried  downstream. 

Computations  have  been  made  simulating  the  nominal  discharge  conditions 
maintained  in  the  Raman  scattering  experiments,  i.e.:  p  -  20  Torr,  E/n  = 

2.5  x  10" -*-5  v-cra2,  i  =  55  mA,  v  =  5  x  10^  cm/sec  and  Tj_  =  300°K.  The  effect 
of  discharge  current  has  been  examined  by  varying  the  discharge  power  den¬ 
sity,  p/v.  Figure  15  shows  the  evolution  of  the  effective  hydrogen  vibra¬ 
tional  temperature  between  the  ground  and  the  first  excited  vibrational 
level  for  several  values  of  the  discharge  power  density  parameter.  The 
nominal  experimental  discharge  current,  i  =  55  mA  corresponds  to  a  power 
loading  of  p/v  =1.7  w/cm^.  The  effective  vibrational  temperature  rises 
rapdily  until  the  end  of  the  discharge  then  decays  at  a  rate  depending  on 
the  power  loading.  The  decay  of  vibrational  excitation  is  mainly  due  to  VT 
decay  by  H  atoms.  Higher  power  loading  leads  to  higher  atom  concentrations 
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Figure  15.  Hydrogen  Discharge/Afterglow  Excitation  Profiles 
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at  the  discharge  exit.  The  H-atom  concentration  at  the  discharge  exit  is 
predicted  to  be  ~  lO-1-?  cm”3  corresponding  to  a  mole  fraction  of  1.6  x  10~3 
percent,  for  P/v  =  1.7  w/cm^. 

The  evolution  of  the  hydrogen  vibrational  distribution  is  shown  in 
Fig.  16  for  a  low  power  loading  parameter.  It  can  be  seen  that  the  effect 
of  the  large  vibrational  anharmonicity  in  hydrogen  mainly  affects  the  vibra¬ 
tional  levels  higher  than  the  fourth.  The  behavior  of  the  lower  levels  is 
overned  mainly  by  BT  energy  transfer  through  atomic  collisions. 

Calculations  have  been  made  down  to  quite  low  discharge  power  loading 
densities  in  order  to  simulate  the  possible  effect  of  large  discharge  ineffi¬ 
ciencies  described  in  the  section  on  discharge  heat  balance  measurements. 

It  should  be  noted  for  P/V  =  0.5  w/cm3,  typical  of  inefficient  discharge 
couplings  measured,  the  first  vibrational  level  temperature  rises  to  1390°K 
at  the  end  of  the  discharge  and  falls  to  3360°K  at  30  cm.  These  values  are 
in  general  agreement  with  the  experimental  results. 
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